Strong coupling and long-range collective interactions in optomechanical arrays 
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We study cavity optomechanics with an array of equidistant scatterers in a Fabry-Perot resonator. 
Our treatment departs from the traditional reflective optomechanics regime with single scatterers 
and introduces a novel transmissive regime that allows better coupling of light to the collective 
motional dynamics of the array. We show that the array can behave as a single, highly transmissive, 
'superscatterer' possessing a specific collective motional mode which is strongly coupled to the cavity 
field with an optomechanical coupling strength that scales as N 3 ^ 2 with the number of scatterers N 
and does not saturate with the reflectivity of the scatterers, as opposed to the reflective case. For 
the case of small polarisability, e.g., many atoms in an optical lattice, our treatment also uncovers 
a new regime where the same jV 3//2 scaling is expected. 



I. INTRODUCTION AND MOTIVATION 

Optomechanics investigates the interaction of light 
with mechanical oscillators [1, 2] ranging from micro- 
scopic ensembles of atoms [3] to macroscopic resonators 
such as nanomembranes, micromirrors, etc. [1, 2, 4-7] 
Besides having implications in fundamental quantum me- 
chanics [2], exploiting motion at the quantum level has 
important applications in metrology [1] and quantum in- 
formation processing [8]. However, despite recent spec- 
tacular progress [1, 2, 9], single-photon coupling is typi- 
cally very weak for massive resonators, necessitating the 
use of many photons to amplify the interaction [6]. We 
introduce the idea that the collective motion of an array 
of scatterers in a Fabry-Perot resonator can couple to the 
cavity field with an interaction strength that scales super- 
linearly with the number of scatterers and does not satu- 
rate with their reflectivity. This new regime should real- 
istically allow for achieving strong single-photon optome- 
chanical coupling [2, 6, 10] with massive resonators and 
realising quantum optomechanical interfaces [7, 8, 11]. 
It also sheds light into the long-ranged cooperative op- 
tomechanical coupling of atoms or membranes in optical 
cavities [3, 12, 13]. 

Cooling a mechanical oscillator to its ground state [9] 
is but one of a series of achievements that demonstrate 
the power of coupling light to moving scatterers. (In this 
article we use terms relevant to the optical domain, but 
it is to be understood that our conclusions hold equiva- 
lently for, e.g., electromechanical systems.) Among the 
various approaches followed to couple mechanical oscilla- 
tors with optical resonators, a successful one involves po- 
sitioning reflecting objects — atoms [3] or dielectric mem- 
branes [4, 5] — inside an optical cavity. In the former 
systems the coupling constant scales as TV 1 / 2 with the 
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number, N, of atoms [3]. In the latter, it saturates to a 
fundamental limit g as the reflectivity of the membrane 
approaches unity [4, 5]. For a highly reflective membrane 
placed close to the centre of a Fabry-Perot resonator of 
length L and resonance frequency uiq, we can calculate 
the single-photon coupling strength as the shift in cav- 
ity frequency when the mirror moves through a distance 
equal to the spread, xq, of its zero-point fluctuations: 
g = 2u!qXo/L. Several approaches are being followed 
to reach the strong optomechanical coupling regime in 
membrane-based systems [4, 5, 14]. 

Meanwhile, interactions between distant elements in ar- 
rays of massive scatterers are believed to be strongly sup- 
pressed [13], whereas atomic systems are said to demon- 
strate infinitely long-ranged interactions [3]. These two 
systems are two limiting cases of a more general model 
we describe below. We identify a remarkable 'trans- 
missive' configuration for the scatterers, regardless of 
whether these are atoms or micromirrors, where the cou- 
pling strength (i) scales as iV 3 / 2 and (ii) does not satu- 
rate as the reflectivity of the elements approaches unity. 
This novel regime allows multi-element opto- or elec- 
tromechanical systems to reach single-photon coupling 
strengths orders of magnitude larger than is currently 
possible. Concomitantly, we show that in this configu- 
ration (iii) the resonator field couples to a 'sinusoidal' 
mechanical mode supporting inter-element interactions 
that are as long-ranged as the array itself. 

Let us again consider a lossless membrane, of thick- 
ness much smaller than a wavelength, placed inside a 
Fabry-Perot resonator. This time, we suppose that 
the membrane has an amplitude reflectivity r, which 
we parametrise in terms of the polarisability £ 

r yj\ — \r\ 2 . The single-photon optomechanical cou- 
pling strength is now g = g\r\, which is maximised to 
g for large i.e., in the reflective regime \r\ — > 1. In 
order to illustrate the emergence of collective optome- 
chanics, we consider now two identical membranes placed 
symmetrically in the resonator at a distance d from each 
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FIG. 1. (Color online) Schematic of our model, motional 
modes, and working points. (a), System considered: N 
equidistant scattering elements positioned in the field of a 
Fabry-Perot resonator, (b), Two examples of collective mo- 
tional modes: the centre-of-mass mode, and a representative 
of a class of 'breathing' modes, (c), Free-space reflectivity of 
a 5-element array as a function of the inter-element spacing d 
(in units of A). The reflectivity of each individual element is 
17.6% (( = —0.5), illustrated by the dotted green line. 



other, in the spirit of Fig. 1(a). As we justify below, 
the effective polarisability of the two-clement system is 
found to be of the form \ — 2£[cos(fccf) — (sin(kd)] for 
light having wavenumber k. This effective polarisability, 
and thereby the reflectivity, vanishes when d is chosen 
such that kd = tan~ 1 (l/C) mod tt. Assuming this trans- 
missive condition, one can linearise the cavity resonance 
condition for a small variation 5d of the mirror spacing. 
This readily gives an optomechanical coupling strength 



9o = 



5u> 



V2t 



1- r 



(1) 



where x' = xq/V^ is the extent of the zero-point mo- 
tion for this breathing mode. It is evident that g' scales 
more favourably with |r| than go. One can interpret 
this result by noting that as the reflectivity of the in- 
dividual elements is increased, the constructive interfer- 
ence that is responsible for making the array transmissive 
also strongly enhances the dispersive response of the cav- 
ity around this working point. In a symmetric situation 
the displacement of a mirror in one direction will cause 
the field to adjust so that the other mirror moves in the 
opposite direction, thereby balancing the power imping- 
ing on the two mirrors. In this simple two-element case, 
the radiation-pressure force thus couples naturally to a 
breathing mode [Fig. 1(b)]. 

Optical ' supers catterers'. — To treat the general case of 
an array of N equally-spaced elements in free space we 
make use of the transfer matrix formalism [15] for one- 
dimensional systems of polarisable scatterers, and derive 
the response of the system to a propagating light field. As 
is well-known from the theory of dielectric mirrors, the re- 
flectivity of the ensemble can be tuned to have markedly 
different behaviors at a given frequency [Fig. 1(c)]. An 
array of N equally-spaced identical elements, each of po- 
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FIG. 2. (Color online) Transmission through a cavity with 
5 elements. The sinusoidal mode, a type of breathing mode 
[Fig. 1(b)], is defined in the text. Note the heightened sensi- 
tivity of the resonance at the centre of the figure. (( — —0.5, 
L — 6.7cm, d — d+, bare cavity finesse w 3 x 10 4 .) 



larisability £, can be described through the transfer ma- 
trix 



M N = 



(1 + ix)e v ix 



which relates left- and right-propagating fields on either 
side of the array. For real £ it can be shown that TV 
lossless scatterers behave as a collective 1 super scatter er' 
having effective polarisability (see Appendix) 



X = C sin [^cos 1 {a)]/\fl 



(2) 



with a = cos(kd) — (sm(kd), together with a phase shift 
/j. The ensemble is attains its largest reflectivity for kd = 

kd = -tan-HC), x = Xo = - i sin[7Vcos- 1 (Vl + C 2 )], 
and becomes fully transmissive (x = 0) for kd = kd± = 
-tan-^C) ± cos" 1 [(l + C 2 r 1/2 cos(7r/A0]. These are 
the two working points illustrated in Fig. 1(c). 

The phase shift fi, which is the phase accrued on re- 
flection from the stack, can be expressed in terms of the 
Chebyshev polynomials of the second kind, U n (x) (see 
Appendix) : 

Au l-i(U N -i(a) 



(l-i()U N ^(a)-e* kd U N _ 2 (a) 



From this point onwards, the array can be treated as a 
single scatterer, keeping in mind the dependence of x on 
the inter-element spacing. 



II. OPTOMECHANICAL COUPLING 
STRENGTH OF THE ENSEMBLE 

When placed inside a cavity, cf. Fig. 1(a), this array 
of scatterers modifies the resonance condition, such that 
the resonances of the system are given by the solutions 
to (see Appendix) 



1 + 



iXCOs{2kx) ± y 1 + x 2 sin (2kx) 



(3) 



These resonances show up when plotting transmission 
spectra of the system, e.g., as the bright curves in Fig. 2 
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FIG. 3. (Color online) Coupling strengths for multi-element 
arrays. All these curves are scaled to g « 2-zr x 15 Hz 
(dotted green line). Top: Scaling, with the number of ele- 
ments N, of the normalised coupling strength for the sinu- 
soidal (oc TV 3 / 2 ) and centre-of-mass (cc TV -1 / 2 ) modes, as 
illustrated by the dotted curves. (£ = —0.5, L = 6.7 cm, 
d — d+ + 20 A.) Bottom: Comparing coupling strengths for the 
sinusoidal (solid lines) and centre-of-mass (dashed) modes, 
with C = -0.5 (17.6% reflectivity; blue) and C = -12.9 
(99.4%; red). (d = d+.) 



in the case of the sinusoidal mode. Generically, one ob- 
tains the linear optomechanical coupling strength by lin- 
earising this equation about one of its solutions. For 
a centre-of-mass motion [cf. Fig. 1(b)] in the reflective 
regime, d = do, we thus obtain g com = g^H/N, where 
7Z — Xo/(l + Xo) * s the maximal power reflectivity of the 
ensemble. As TV or £ increase, 1Z saturates to 1 and gcom 
scales as TV -1 / 2 . This scaling can be explained simply 
by noting that the motional mass tun of TV elements is 
TV times that of a single one; the single-photon coupling 
strength, which is proportional to 1/^/mjv, therefore de- 
creases with TV. 

In the transmissive regime, d — d + , the intensity profile 
peaks at the centre of the array [Fig. 1(a)]. The optome- 
chanical coupling strength gj for the j th membrane is 
strongest where the difference in amplitudes across the 
membrane is greatest, j ps (TV + 2)/4 or (37V + 2)/4, re- 
sulting in gj oc sin[7r(2j — 1)/7V]. The cavity field there- 
fore couples to the sinusoidal mode, whose profile varies 
sinusoidally across the array, with a collective coupling 
strength (TV > 2) (see Appendix) 



fN 3lC|csc(£)^sin 2 (£)+C 2 + |CI 
2 l + 2/Vf|C|csc2(^) v /sin 2 (^)+C 
or, for large TV, 

^g( 2 N 3 / 2 V2 
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FIG. 4. (Color online) Linewidth-narrowing effect. The effect 
becomes significant as the reflectivity of the individual mem- 
branes increases, and is more pronounced for widely-spaced 
membranes. Inset: Representative transmission profiles for 
N = 10 for the parameters in the main figure. (L = 6.7 cm.) 



the last expression is valid for L/d ^> 2( 2 N 3 /tt 2 . Op- 
timising over TV for arbitrary L/d, we obtain g°? = 
\g^J L/d\C\- This favourable scaling with both TV and 
|r|, as shown in Fig. 3, is a significant improvement over 
the state of the art. Equation (4) is valid both for ar- 
rays of membranes (|£| > 1) and of atoms <C 1). 
For typical experiments with cold atoms in cavities [3], 
N\C\ <C 1, i.e., the atoms do not significantly modify 
the mode structure of the cavity resonance. In this case, 
f/ S i n reduces to the usual TV 1 / 2 scaling that arises from 
the independent coupling of well-localised atoms inter- 
acting with an unmodified cavity field [3]. However, our 
approach allows us to treat the case of a confined, opti- 
cally dense [16], and well-localised atomic ensemble hav- 
ing N\C\ > 1 [17], which would exhibit a superlinear 
TV 3 / 2 scaling, signalling the breakdown of independent 
coupling. 

An interesting effect arises in the regime where gsin 
saturates and eventually starts decreasing as a func- 
tion of TV; the scatterers then act to narrow the cav- 
ity resonance substantially. This has possible applica- 
tions in hybrid systems, along the same lines as those 
of electromagnetically- induced transparency in Ref. [11], 
and arises from an effective lengthening of the cavity, 
due to the presence of the array, to a length L c g = 
L+ ^dC, 2 N 3 . From the fact that the cavity finesse in the 
transmissive regime is fixed by the end mirrors, it follows 
that the linewidth of the cavity is K c g oc 1/L c g, cf. Fig. 4. 
When using low-finesse cavities and low mechanical oscil- 
lation frequencies, this effect could thus be used to place 
the system well within the sideband-resolved regime. 



III. NUMERICAL EXAMPLE 

The power of this approach to optomechanics is best 
seen through a numerical illustration. If we take commer- 
cial silicon nitride membranes [4] with a power reflectiv- 
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FIG. 5. (Color online) Matrix of inter-element forces. The 
long range of the interaction between elements of the array 
can be illustrated by calculating the force experienced by any 
element (horizontal axis) due to the motion of any other ele- 
ment (vertical axis); for a symmetric system such as ours this 
matrix is also symmetric. The strongest interactions corre- 
spond to the darkest blue (in-phase) or orange (out-of-phase) 
regions. 

ity of 17.6% (C = —0.5) and xo — 1.8 fm, and a cavity 
with L = 6.7 cm and ojq equivalent to a wavelength of 
1064 nm, we can estimate g com ps 2n x (12.8 x N^ 1 / 2 Hz) 
for N > 3. For the sinusoidal mode, and with the same 
parameters, g S i n rj 2tt x (1.3 x N 3 / 2 Hz) for large N; 
an improvement by over an order of magnitude when 
N = 10. A transparent ensemble potentially provides 
a much stronger optomechanical coupling than a reflec- 
tive one; indeed (fen/Scom cx iV 2 . Let us now consider 
highly-reflective membranes [14] having 99.4% reflectiv- 
ity (£ = -12.9), x = 2.7fm, and uj m = 27rx211kHz. For 
L = 0.25 cm and d = d +1 one obtains g CO m ~ 27r x 600 Hz 
and g S i n R3 2-7T x 270kHz > ui m (N = 5); strong coupling 
between a single photon and a single phonon is already 
within reach with only very few elements. 

Long-range collective interactions. — The collective na- 
ture of the interaction that is responsible for these large 
coupling strengths also gives rise to an effective 'non- 
local' interaction between the scatterers, where the mo- 
tion of any particular element influences greatly elements 
further away, and not just its nearest neighbours. The 
cavity field mediates an effective interaction between 
pairs of elements (indices I and j, say) through a cou- 
pling constant proportional to gigj, illustrated as a 'force 
matrix' in Fig. 5. This coupling is macroscopically long- 
ranged, and our system therefore provides an ideal means 
to gain insight into collective optomechanics phenom- 
ena [13]. By contrast, in the reflective regime the light 
does not permeate through the ensemble, and the force 
matrix would instead be almost diagonal. 

IV. LIMITATIONS 

Our model makes three key assumptions about the 
physical make-up of the membrane ensemble, which we 
will now discuss in brief, (i) The membranes are as- 



sumed to be thin on the scale of a wavelength. Since 
silicon nitride membranes have thicknesses of the order 
of 50 nm <C A ~ 1 itm [4] , this assumption is not expected 
to be a limiting factor. 

(ii) Positioning errors are also a concern. Piezoelectric ac- 
tuators can be used to position the membranes with sub- 
nm accuracy. Nevertheless, let us estimate the position- 
ing requirements in a fairly demanding setup [N = 10, 
£ = —0.5). The reflectivity of the ensemble reaches 10% 
when the spacing between each pair of membranes is 
off by ~ 10~ 3 Awl nm. This estimate is conservative, 
since it assumes that the spacing is uniformly incorrect 
throughout the entire ensemble. If the error in the posi- 
tion of any element is independent of that of any other 
such that the effects add up 'incoherently,' the expected 
tolerance required is not beyond current technology. 

(iii) The membranes were also assumed to be non- 
absorbing; this is an excellent approximation to a sin- 
gle membrane with an imaginary part of the refractive 
index being < 10 -6 -10 -5 [5, 18]. In the transmissive 
case, the optical field permeates throughout the entire 
ensemble; it is precisely this fact that gives rise to the 
collective dynamics we described above. The power ab- 
sorbed by the ensemble therefore increases (oc iV~ 2 ) with 
N, and even moderate absorption can seem prohibitive 
for large ensembles of membranes. This is mitigated by 
the large coupling strengths obtained, which allow much 
smaller photon numbers to be used: g 2 ln oc TV 3 increases 
faster than the power absorbed as N increases. We note 
further that at large input powers it might be possible 
to exploit the photothermal force to further enhance the 
collective optomechanical interaction [14, 19, 20]. 



V. CONCLUSIONS 



We have shown that by operating in the transmis- 
sive regime of an ensemble of polarisable scatterers, 
one couples to a collective mechanical mode strongly 
enough to bring single-photon optomechanics [10] and 
strong-coupling effects [6] within reach. The system also 
gives rise to inter-element interactions that span the en- 
tire array. Our ideas apply generically across a wide 
range of systems; any system that can be modelled as 
a one-dimensional chain of scatterers (e.g., membranes, 
atoms [16], optomechanical crystals [9], or dielectric mi- 
crospheres [21]) is amenable to a similar analysis and 
shows the same rich physics. Similar methods would al- 
low the extension of these ideas to more complicated sys- 
tems where the polarisability is a function of frequency 
or of position along the array, or systems involving the 
interaction of arrays of refractive elements with multiple 
optical modes. 
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Appendix A: Transfer matrix for iV-membrane stack 

Let us consider N equally-spaced identical non- 
absorbing membranes, each of which has a thickness 
much smaller than an optical wavelength. The spacing d 
between these membranes determines the overall optical 
properties of the ensemble. We shall treat the system as 
a strictly one-dimensional system, and we shall use the 
transfer matrix formalism [15, 22]. Our starting point is 
the matrix that links the fields interacting with a single 
membrane of polarisability ( (( e R for a lossless mem- 
brane), 



M m (C) = 



l + »C K 
-K i-*C 



(Al) 



and the matrix that describes propagation of a 
monochromatic beam of wavenumber k over a distance d 
through free space, 



M p (d) EE 



/A''/ 









— ikd 



(A2) 



both of which have unit determinant. We wish to evalu- 
ate a product of the form 



M m (0 ■ M p (d) ■ M m (0 ■ ■ ■ M m (() , 
with N copies of M m (Q. First, we note that 

M p (d/2) ■ M m (C) • M p (d) ■ ■ ■ M m (C) • M p (d/2) 
= [M p (d/2) ■ M m (C) • M p (d/2)] N ee M N . 



(A3) 



where the second line defines the matrix M: 
M 



(l + iQe lkd i( 

-i( (1 - i()e- lkd 



(A4) 



(A5) 



We can easily see that detM = 1, whereby it can be 
shown [23] that for real £ we can write 



M 



rN 



(i + i X y {kd+ti) 

-iX (1 



-i(kd-\-fj,) 



(A6) 



where x = C^v-i( a )i with U n {x) being the n th Cheby- 
shev polynomial of the second kind, a — cos(kd) — 
(sin(kd), and 



1 - i(U N -i(a) 



(1 - iC)C^-i(o) - e ikd U N _ 2 (a) 



(A7) 



Upon removing the padding of d/2 from either side, we 
obtain the matrix that describes the TV-membrane en- 
semble: 

M P W (2k)} ■ M m ( X ) ■ M p [^/ (2k)) . (A8) 

Thus, as stated in the main text, N lossless membranes 
behave as a collective 'supermembrane' of polarisability 
X along with a phase shift /i/2 on either side of the stack. 



Appendix B: Resonances of the compound system 

Following the main text, let us now place our mem- 
brane stack inside a perfect Fabry-Perot cavity. Finding 
the resonances for the compound optical system is seem- 
ingly trivial; we need only solve the relation 







-iO 



1 



-%X ix 
i x 1 - ix 
i<f> o 







(Bl) 



with 6 ee k(L/2 + x) + /i/2 and <p = k(L/2 - x) + pt/2. 
Here x is the displacement of the ensemble with respect 
to the centre of the cavity. We thus obtain 



ikL 



1 + iX 



iX cos(2/ca;) ± y 1 + x 2 sm (2kx) 



(B2) 



However, we immediately see that this equation is tran- 
scendental in k, and therefore cannot be solved analyti- 
cally; this equation is easily solvable for L, and we there- 
fore know the resonant length of our cavity. It is now a 
legitimate question to ask: 'If d (or x) shifts by a small 
amount, how much will the resonant frequency of this 
cavity shift?' This question is, of course, easily answered 
by expanding Eq. (B2) in small increments about its so- 
lution. Assuming a dominantly linear effect, we replace 
k — » ko + Sk, x — > x + Sx, x — > X + an d M — > M + 
in Eq. (B2). Around resonance, the result simplifies to 



L8k + <5/i = 



-1 ± cos(2k x) J y 1 + x 2 sm 2 (2k x) 

x 6 X /{l + X 2 ) 
2x sm(2k x) j \Jl + x 2 s\n 2 (2k x) 

x (xSk + k Sx) . (B3) 



() 



Appendix C: Centre-of-mass coupling strength 

Let us start from Eq. (B3). For the centre-of-mass 
motion, dfj, — d\ = 0, and we assume that \L/x\ is very 
large, such that we can write 



LSk = T 



2xsm(2k x) j \J I + x 2 sin 2 (2fc a;) 



k Sx . 



( C1 ) 

The right-hand-side of this equation is maximised when 
sin(2fcoa;) = qpl, whereby 



L5k = 2k (-x/Vl + X 2 )6x- 



(C2) 



This is, in absolute value, a monotonically-increasing 
function of \x\ and is therefore maximised when x at- 
tains its largest value, ( Un-i + C 2 )> which leads to 
g CO m as defined in the main text. 



Appendix D: Coupling to each individual membrane 

The matrix Mn representing the ensemble can be writ- 
ten, for 1 < j < N, 



e -Vi/ 2 J [ -ixi 1 - ixi 

e —i(ni/2+u+kxj) 

1 + iC K 1 \ e i ^/' 2+v - kx ^ 

— iC, 1 — i( e -i(ti2/2+i/-kxj) 

1+% % 1 |"e V2/2 
-%Xi 1-«X2 e~ V2 / 2 



, (Dl) 



where /Ui and xi describe the ensemble formed by the 
ni = j — 1 membranes to the 'left' of the j th , and (X2 and 
X2 the one formed by the n 2 = N — j membranes to its 
'right'. The displacement of the j th element is denoted 
Xj ; all other membranes arc in their equilibrium position. 
In the transmissive regime, to lowest order in kxj in each 
entry, the matrix product above can be written, with the 
above choice for 



axi 



j3* Xj 



fiXj 

+ a*Xj 



(D2) 



where a and /? are increments of first order in the relevant 
displacement [note that the (off-)diagonal terms are com- 
plex conjugates of each other; this is different to the case 
where absorption is nonzero] . When this matrix is substi- 
tuted into the equation for the resonance condition, the 
terms involving Re{e _IA1 a} and Re{/3} drop out entirely 
for a symmetric system, such that it suffices to consider 
only the imaginary part of the increment. Let us reiter- 
ate that this happens only because the off-diagonal terms 
are complex conjugates of each other; were absorption to 



be nonzero, this would no longer be the case. Eq. (B3) 
now simplifies to 



dk _ Im{P + e-^a} 
dx~~ Z+2dE 



(D3) 



with v = kd, 
a = 2ik( 

= 2ik( 2 



and 



P = 2kQ 



e Vl (l + %)X2-e V2 Xi(l + %) 

(l + e)U 2 ni -i{a)Un 2 -i{a) 
(1 - iOU nx -i(a) - e iv U ni - 2 (a) 

(l-iC)C/„ 2 -i(a)-e^[/„ 2 _ 2 (a) 



XiX2-(l + »Xi)(l-%)e l(AI1 ^ 2) 



(D4) 



2fcCK 2 C4 1 -i(a)[/„ 2 -i(a) 



(l-i()Un 2 -i{a)-e™U n2 - 2 (a) 
X (1 - iC)U ni -i(a) - e iv U ni - 2 {a) 



(D5) 



These two expressions simplify considerably to yield 



Im{^ + e~ v a} = 2fcCcsc(-^ 



n) 



^n-j + c 2 -c 



sin 27r 



N 



(D6) 



and therefore 



9j 



-2uj c x 



Ccsc(f) Jsin 2 (^)+C 2 -C 



L - 2Nd( esc 2 (%) ^/sin 2 (^)+C 2 



x sin 2-7T 



N 



(D7) 



As discussed in the main text, the coupling of the collec- 
tive motion of the membranes to the cavity field is gov- 



erned by the constant \jY^j=i 9j> such that for N > 2 



N ffCcsc(f, 



sin 



1( tT 



+ c 2 -c 



!-2iVfCcsc 2 (^ 



+ c 2 



because of the relation 



N 



^sin 2 (2^ 

\ 3 = 1 



N 



V2 for N = 2 



for TV > 2 



(D8) 



(D9) 



this is equal to g s i n as defined in the main text in the 
appropriate limits. 
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